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^ : ABSTRACT 
O 

o . 

^ ' The recent detection of variable infrared emission from Sagittarius A*, com- 

■ bined with its previously observed flare activity in X-rays, provides compelling 
evidence that at least a portion of this object's emission is produced by non- 

Q> ' thermal electrons. We show here that acceleration of electrons by plasma wave 

turbulence in hot gases near the black hole's event horizon can account both for 
Sagittarius A*'s mm and shorter wavelengths emission in the quiescent state, 
and for the infrared and X-ray flares, induced either via an enhancement of 

■ the mass accretion rate onto the black hole or by a reorganization of the mag- 

ff^ ■ netic field coupled to the accretion gas. The acceleration model proposed here 

O 

O ■ multi-wavelength observations. We further suggest that the diffusion of high en- 

^ ■ ergy electrons away from the acceleration site toward larger radii might be able 

to account for the observed characteristics of Sagittarius A*'s emission at cm and 



> 



2 ■ longer wavelengths. 



Subject headings: acceleration of particles — black hole physics — Galaxy: center 
— plasmas — turbulence 



1. Introduction 

The observation of stellar motions within light-days of Sagittarius A*, a compact radio 
source at the Galactic Center (Balick & Brown 1974), has provided compelling evidence that 
this source is the radiative manifestation of a ~ four million solar mass black hole (Schodel 
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et al. 2002; Ghez et al. 2003). The recently detected infrared emission and flare activity 
have provided an additional evidence that this source is powered by a hot gas accreting onto 
the black hole (Baganoff et al. 2001; Goldwurm et al. 2003; Baganoff 2003; Porquet et 
al. 2003; Zhao et al. 2004; Ghez et al. 2003). The quasi-periodic near-infrared variability 
may be an indication that the gas flared up before spiraling into the black hole (Genzel et 
al. 2003). It is now generally agreed that the radio and infrared emission, and the flares, 
are hkely produced by nonthermal high energy electrons (Liu & Meha 2001; Genzel et al. 
2003; see also Mahadevan 1998). However, the exact nature of the mechanism responsible 
for the acceleration of the electrons has not been addressed. This has given rise to diverse 
interpretations of the observations with assumed spectra of the accelerated electrons (Markoff 
et al. 2001; Liu & Meha 2002a; Nayakshin et al. 2003; Yuan et al. 2003). 

In this letter, we show that the mechanism producing high-energy particles in solar 
flares works equally well in hot plasmas near the black hole. The solar flare model is based 
on a second order Fermi acceleration process or a stochastic acceleration (SA) of particles 
by interacting resonantly with plasma waves or turbulence (PWT) generated via an MHD 
dissipation process (see e.g. Miller & Ramaty 1987; Hamilton & Petrosian 1992; Petrosian 
& Liu 2004). In Sagittarius A*, nonthermal particles can be produced by the turbulence 
expected to be induced by the magneto-rotational instability in the accretion torus (Balbus 
& Hawley 1991; Melia, Liu & Coker 2001). The radiation emitted at the acceleration site 
can explain the quiescent state mm and shorter wavelength observations. Solar flares are 
energized by the process of magnetic reconnection during the dynamical evolution of the 
coronal magnetic fleld. Similar processes in Sagittarius A* can release energy in a small 
region and produce what we call a local event. A global flare can be induced by an MHD 
fluctuation in the accretion torus or an enhancement of the accretion rate onto the black 
hole. The emission spectra from these two energization mechanisms are quite different, and 
can explain the distinct flare behaviors. 

In § 2, we outline the theory of SA. Its application to Sagittarius A* is presented in § 3. 
The main results of this letter are summarized in § 4, where we also discuss consequences of 
energetic electrons escaping the acceleration site. If such electrons can diffuse toward larger 
radii, they may account for Sagittarius A*'s emission at cm and longer wavelengths and has 
the potential to explain the observed linear and circular polarization characteristics, and 
variabilities of this source. 



-3- 



2. Stochastic Electron Acceleration 

In the theory of SA (Petrosian & Liu 2004), particles are accelerated from a background 
plasma to higher energies by interacting resonantly with PWT. The particle distribution and 
the consequent emission spectrum are determined by the acceleration (-Dee) and scattering 
(D^fj^) rates due to the PWT, the energy loss rate and the spatial diffusion time Tesc- 
These rates and time are determined by the magnetic field B, density n and particle distri- 
bution of the background plasma, the spectrum and intensity of the turbulence, and the size 
of the acceleration site R. 

The wave modes in a plasma are described by the dispersion relation cu — a;(k), which 
depends primarily on the plasma parameter a — ujpe/^e — (47rmen)^/^c/S , the ratio of the 
electron plasma frequency ujpe — {4:TTne^ /meY^"^ to the nonrelativistic electron gyrofrequency 
Q,e — eB/nieC, where rric, c and e are the electron mass, the speed of light, and the elemental 
charge unit, respectively. Electrons with Lorentz factor 7, velocity v, and pitch angle cosine 
H couple strongly with the wave satisfying the resonance condition: uj = k\\v/j + j^^ , where uj 
and k\\ are the wave frequency and wave vector parallel to the large scale magnetic field (only 
the first harmonic of the gyrofrequency is considered here). Electrons with higher energies 
resonate with waves with smaller wave numbers (k) corresponding to larger spatial scales. 

In the following discussion, we adopt a turbulence spectrum characterized by a broken 
power law with a lower wavenumber cutoff at /Cmin, corresponding to the scale where the tur- 
bulence is generated (< -R), and a break wavenumber k^ax above which the wave dissipation is 
fast. We will study the SA by a PWT propagating parallel to the large scale magnetic field. 
The characteristic interaction time Tp is given by — (7r/2)Jle(87r£^tot/-B^)(? — l)^mTn , 
where q > 1 is the turbulence spectral index in the inertial range and £tot is the total 
turbulence energy density. Relativistic electrons interacting with waves in a plasma with 
a > (mc/mp)^/^ ^ 1/43, where nip is the proton mass, will have an isotropic distribution 
because of the short scattering time = l/D^^. In such plasmas electrons can gain energy 
by interacting with the waves and lose energy either via Coulomb collisions with background 
plasma particles (at lower energies) or by radiative processes (at higher energies). The elec- 
trons also diffuse spatially and leave the region on a characteristic time Tesc ~ / 
The spatially integrated electron distribution N{E), as a function of the kinetic energy E, 
satisfies the well known diffusion-convection equation: 

dN d'^ d ■ N 

where Dee and A are obtained from the wave-particle interaction rates and Q is a source 
term. In Sagittarius A*, where the photon energy density is at least one order of magnitude 
lower than the magnetic field energy density, the radiative loss of relativistic electrons is 



-4- 



dominated by synchrotron process, for which 

Tioss = E/E^ = 47rr2cn(lnA//3 + 4/3V/9«') ' 

where In A ~ 20 and the first and second term in the denominator give the Coulomb coUision 
and synchrotron loss rate, respectively. 



3. Application to Sagitteirius A* 

To determine the exact turbulence spectrum, one needs to solve the coupled kinetic 
equations of the waves and particles and their solutions depend on the wave generation, 
cascade, dissipation processes and the distribution of the background particles (see e.g. 
Miller, LaRosa & Moore 1996). The lack of a well established theory for the MHD turbulence 
renders such a treatment unpractical. As an approximation, we set /Cmin = 27r77/it! with 77 > 1, 
where the characteristic variation length of the large scale magnetic field R/rj is related to 
the turbulence generation length scale. To estimate /Cmax, we note that the wave damping is 
efficient when the waves start to resonate with the background particles. From the resonance 
condition, one has fcmax ~ ^c/{c < 7 >), where < 7 > is the mean Lorentz factor of the 
background electrons, which is expected to be ~ 50 near the black hole of Sagittarius A*. 
For /Cinin < k < /Cmax the turbulence most likely has a Kolmogorov spectrum with an index 
q — 5/3 and above /Cmax the index changes to 4.0, a typical value obtained from MHD 
simulations (Hawley et al. 1995; Vestuto et al. 2003). MHD simulations (Hawley et al. 
1995; Hirose et al. 2004) also suggest that the magnetic field energy density is about a few 
percent of the energy density of the background gas, which means that a ~ 1 (note that 
= 0.5nmeC^/[-B^/87r]), and that Stot ^ -B^/Svr. The accretion fiow in Sagittarius A* is 
coUisionless; it is not clear how the particle distribution of the accretion plasma evolves from 
large to small radii, where it enters the acceleration site of strong dissipation and PWT. 
We make the reasonable assumption that the source electrons have a relativistic Maxwellian 
distribution. We therefore have £tot = o,nk]jT with a <C 1 as a model parameter. The exact 
spatial diffusion time (or Tgsc) depends on the structure of the magnetic field (Hirose et al. 
2004) and the scattering rate. We set Tesc = '^ul'^sc + 7"tr, where the transit time Ttr = rjR/v 
determines the escape time when the scattering time Tgc ^ Ttr (Petrosian & Donaghy 1999; 
Petrosian & Liu 2004). 

During big flares, the physical conditions are expected to undergo dramatic changes, 
which may affect the turbulence spectrum and the energy partition among the gas, turbu- 
lence, and magnetic field. These processes are not well understood. For the purpose of 
illustration, we will make the following reasonable assumption to reduce the number of pa- 
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rameters. We set = 2, k^^ = 0.02 Qe/c = 1.2 x 10-^{B/G) cm-\ a = 1 and a = 0.008. 
This leaves only three independent parameters; here we choose R, n and Tp. All other 
quantities can be derived from them. 

To calculate the emission predicted by the model, one also needs to specify the geometry 
of the source. We will assume that the source is uniform and spherically symmetric with 
a radius R, which may or may not be centered on the black hole. (Note that although 
the actual geometry is inferred to be Kcplerian for a global acceleration region — see, e.g., 
Melia et al. 2001a — its differences from the simplified geometry adopted here are expected 
to change the required model parameters only slightly. This is a fair assumption for local 
fluctuations.) For flares induced by an enhanced accretion rate, the density n can increase 
dramatically while R and Tp may not change much at all. For local MHD fluctuations, all 
three parameters can experience signiflcant changes (Hawley et al. 1995). 

In Table 1, we summarize parameters of several models that account for Sagittarius A*'s 
spectrum in its various states (Figure 2). The corresponding time scales and normalized 
electron distributions in the steady state are shown in Figure 1 for Models A (left panel), 
B (middle panel) and C (right panel). The valleys in the acceleration and the peaks in 
the escape time correspond to the spectral break /cmax of the turbulence. The loss time 
Tioss = E/Ej^ peaks at an energy, which only depends on the plasma parameter a (see eq. [2]). 
In the lower panels of the figure, we also show the source thermal electron distribution and 
the fiux / = N/Tcsc of escaping electrons. It is crucial to note that the acceleration time for 
Model B (Model C) is shorter than the duration of X-ray flares with a hard (soft) spectrum. 
We can therefore use the steady state electron distribution to calculate the synchrotron and 
synchrotron self-Comptonization (SSC) photon spectra of the models. These are compared 
with Sagittarius A*'s spectra in its various states in Figure 2, where the data are gathered 
from observations made at different epochs (Falcke et al. 1998; Zhao et al. 2003; Miyazaki 
et al. 1999, 2003; Serabyn et al. 1997; DoweU 2003; Cotera et al. 1999; Stolovy et al. 
2003; Baganoff et al. 2003). 

Model A explains Sagittarius A*'s emission at frequencies above 100 GHz in the quies- 
cent state. The mm to infrared emission is produced via synchrotron processes, whereas the 
optical to gamma-ray emission is produced by SSC (see Melia et al. 2000). Global fluctua- 
tions in the PWT will vary Tp and result in a different spectrum. For example, by increasing 

by 5%, we obtain Model A", which predicts a weak flare with a soft spectrum. The 
flux densities in the infrared and X-ray bands increase by ~ 80%. Such a flare is clearly de- 
tectable in the infrared by VLT and Keck, but not by the existing X-ray instruments due to 
its low flux. If the flare only lasted for a few hours, neither Chandra nor XMM-Newton would 
be able to detect it. Actually, Sagittarius A*'s quiescent state X-ray spectrum is obtained 
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by averaging over a long observation period (Baganoff et al. 2003), which will smooth out 
contribution from such weak flares^. This aspect of the model therefore may explain why the 
occurrence frequency of infrared flares is more than twice that of X-ray flares (Genzel et al. 
2003). 

Model B (B") reproduces X-ray (and probably infrared) flares with hard spectra. 
Compared with Model A (A"), it has a much smaller radius, which suggests that this is a 
localized (as opposed to global) event. The quiescent state emission of Model A presumably 
still exists during this type of flare. If the plasma is in a Keplerian motion (see Melia et al. 
2001), one would expect to detect quasi-periodic variability, as has been seen in the infrared 
(Genzel et al. 2003), and apparently now in X-rays as well (Aschenbach et al. 2004). Note 
that in the energy range 80 < 7 < 2000 the acceleration time is shorter than the escape time 
(middle panel of Figure 1). This means that the acceleration of electrons in this range is more 
efficient than at lower energies, resulting in a harder electron distribution and corresponding 
harder photon spectrum. These spectra are cut off sharply where tioss, T'esc < Ta- Such flares 
are likely produced by a local fluctuation (possibly a reconfiguration or reconnection of the 
magnetic field), which produces strong PWT, trapping electrons within the acceleration site 
more efficiently. This model explains the 10-27-2000 flare observed by Chandra and predicts 
strong infrared emission accompanying the X-ray flare. Model B" is obtained by decreasing 
R hy a, factor of 13 while keeping nR oc (TaTesc)^'^^/Tioss and R/r^ oc {T^d'Ta)^^'^ the same 
as in Model B, so that the spectrum does not change significantly. This model predicts 
weaker infrared and X-ray flares, but the spectra are harder, which may explain the weak 
X-ray flares observed by Chandra (Baganoff 2003), and the recently observed infrared flares 
(Genzel et al. 2003; Ghez et al. 2003). 

In Model C, which explains the 10-3-1002 XMM-Newton flare, the synchrotron loss 
dominates at very low energy (7 ~ 200). The corresponding electron distribution and X-ray 
spectrum are steeper. Given the required high density and large volume for this model, this 
type of event is likely induced by an enhancement of the accretion rate onto the black hole. 
Future multi-wavelength observations will be able to test the high flux of sub-mm and IR 
radiation predicted by this model. 



^ Just before the submission of this paper, wc learnt that coordinated IR and X-ray observations seem to 
have detected a similar flare (GCNEWS Vol.17; http://www.aoc.nrao.edu/ gcnews) 
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4. Summciry and Discussions 

Based on the theory of stochastic particle acceleration, we have built a model to account 
for Sagittarius A*'s emission at mm and shorter wavelengths. The quiescent state emission 
is attributed to electrons accelerated by turbulence in a magnetized accretion torus and 
flares can be produced via two distinct mechanisms. The IR and X-ray flares with harder 
spectrum are hkely induced by a local MHD process, such as magnetic reconnection. Global 
fluctuations generally produce flares with IR and X-ray spectral indexes close to their cor- 
responding value in the quiescent state. The model not only accounts for the varied spectra 
of the flares, but also explains their relative occurrence rate observed at IR and X-ray. 

Radio emission at longer wavelengths cannot be produced within such a small emission 
region (Liu & Meha 2001). This is not surprising, given that the radio emission is variable 
on a longer time scale of tens of hours to one week, suggesting that this radiation is produced 
at larger radii than what we have been considering here (Zhao et al. 2004, 2003; Zhao & 
Goss 1993; Bower et al. 2002). Very interestingly, our model suggests a significant outflow 
of high-energy electrons (Liu & Melia 2002b), which may very well be the particles that 
eventually produce Sagittarius A*'s cm spectrum via synchrotron emission on a spatially 
larger scale. In the quiescent state (Model A), the power carried away by electrons with 
7 > 100 (electrons with lower energy may be trapped by the gravitational potential of the 
black hole and will not diffuse to larger radii) is about 2 x 10^^ ergs s~^, which is more than 
enough to power the observed radio emission, whose luminosity is about IQ-'^"^^ ergs s"-*". 
(Note that in this model, the mass accretion rate is ~ 10^^~^^g s~^, for which a few percent 
of the dissipated gravitational energy is carried away by the outwardly diffusing electrons; 
see the caption of Table 1.) For the XMM-Newton flare on October 3 2002 (Model C), the 
total energy carried away by the escaping high-energy electrons is about 7 x 10*^° ergs which 
for a few days can sustain the radio flare observed 13 hours after the X-ray event (Zhao et 
al. 2004). This picture also explains the weak correlation between the Chandra X-ray flares 
and the radio emission from Sagittarius A* (Baganoff 2003), because the flux of high energy 
electrons produced during these flares is much smaller than that for the XMM-Newton flare; 
no strong enhancement of radio emission is expected following a Chandra-type of X-ray flare. 

One of the intriguing properties of Sagittarius A*'s radio emission is that it is circularly 
polarized below 100 GHz, even though no linear polarization has been observed there (Bower 
et al. 2002). However, in the mm and sub- mm range only strong linear polarization is 
observed (Aitken et al. 2000; Bower et al. 2003). The mm/sub-mm polarization is 
likely associated with the structure of the magnetic field at small radii (Agol 2000; Melia 
et al. 2000; Bromley et al. 2001), which is consistent with the spectral formation at 
these wavelengths discussed in this paper. The observed circular polarization may be due 
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to anisotropy of the escaping electrons as they are transported along magnetic field lines 
toward larger radii under the influence of synchrotron losses (McTiernan & Petrosian 1990). 
Relativistic electrons beamed along magnetic field lines produces synchrotron emission with 
significant circular polarization while the degree of linear polarization can be suppressed by 
irregularities in the source magnetic field (Epstein 1973; Epstein & Petrosian 1973). Work 
to demonstrate this feature self-consistently is in progress, and the results will be reported 
elsewhere. The model we have presented here shows promise in being able to account not only 
for the spectral characteristics of Sagittarius A*, but also for its polarization characteristics 
and time variation properties. 

Finally, wc emphasize that although we suggested that rj, /Cmax/^e; ct, and a may not 
change significantly over time, and have fixed them in all the models discussed above, given 
the dramatic changes of the physical conditions during a flare state, these parameters may 
also vary slightly, giving rise to different emission spectra than those presented here. Indeed, 
if the infrared upper limits reported by Hornstein et al. (2002) arc confirmed, one then 
needs to decrease the high frequency cutoff of the synchrotron emission by adjusting these 
parameters. This will be investigated in a more comprehensive study of particle acceleration 
in Sagittarius A*. 
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Fig. 1. — Top Panels: The time scale for acceleration = E/A (solid lines), escape 
(dashed lines) and loss rioss = E/Ei^ (dotted lines) for Models A (left panel), B (middle 
panel) and C (right panel). The model parameters are given in Table 1. The loss term 
is dominated by Coulomb collisions below the peak and by synchrotron losses above it. 
Bottom Panels: The normalized electron distributions N (thick lines) at the acceleration 
site and the corresponding escape fluxes / = N/T^sc (thin lines). The long dashed lines show 
the source Maxwellian electron distribution. See text for details. 
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Fig. 2. — Model fits to the broadband spectrum of Sagittarius A*. The data are gathered 
from observations made at different epochs (the instruments are indicated in the figure by 
the numbers). In the radio and infrared bands, the upper set of data show Sagittarius A*'s 
peak fiux densities during the brightest fiares observed so far, while at the same frequencies, 
the lower data points correspond to its emission in the quiescent state. The middle set of 
data in the radio are the averaged flux densities. The upper butterfly in the X-ray band 
gives the peak flux density for the flare observed by XMM- Newton on October 3 2002. The 
middle one gives the peak flux density for the 10-27-2000 flare observed by Chandra. The 
lower one corresponds to the averaged X-ray emission in the quiescent state. The upper 
limits are also for the quiescent state emission. The spectrum produced by each of the flve 
models in Table 1 is indicated by a continuous curve. See text and the caption of Table 1 
for details of the models. 
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Table 1: Models and Parameters. 



Models 


R(rs) 




n (10'^ cm ^) 


B (Gauss) 




Q (lO^Vi) 


A (A") 


2.5 


0.74 (0.78) 


0.76 


8.8 


26.0 (27.3) 


3.2 (3.0) 


B (B") 


0.22 (0.017)t 


26.7 (345) 


19.0 (240) 


44 (160) 


34.2 (20.5) 


0.11 (0.014) 


C 


1.9 


1.78 


18.0 


43 


10.5 


49 



'^The radius R is the radius of the assumed spherical emitting region. In some models, R < rs, meaning that 
the fluctuation producing the flare is not global, but rather is localized near the black hole. See text for 
details. 

Note. — A sampling of models for Sagittarius A* and their corresponding parameters. Note that the first 
three parameters are the primary parameters. All other quantities are derived from them. Q is the steady 
state injection rate. For Models A (A") and C, the mass accretion rate onto the black hole can be estimated 
as M ~ QrUp because high energy electrons (7 > 100), which may diffuse to larger radii, account for less 
than 10% of Q. Sec text for details. Models A, B and C correspond to the left, middle and right panels 
of Figure 1 and the thick solid, dotted and dashed curves in Figure 2, which fit Sagittarius A*'s emission 
in the quiescent state, at the peak of the 10-27-2000 Chandra flare and the 10-3-2002 XMM-Newton flare, 
respectively. The parameters for model A" and B" are indicated in the parentheses whenever different. The 
corresponding emission spectra are indicated, respectively, by the thin solid and dotted curves in Figure 2, 
which explain the weaker infrared and X-ray flares. 



